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TECHNICAL NOTE 3568 


ANALYSIS OF BEHAVIOR OF SIMPLY SUPPORTED FLAT 
PLATES COMPRESSED BEYOND THE BUCKLING 
LOAD INTO THE PLASTIC RANGE 
By J. Mayers and Bernard BudlansEy 

SUMMARY 


An analysis is presented of the postbuckling behavior of a simply 
supported, square flat plate with straight edges compressed beyond the 
bucEling load into the plastic range. The method of analysis involves 
the application of a variational principle of the deformation theory 
of plasticity in conjunction with ccmiputations carried out on a hi^- 
speed calculating machine. Numerical results are obtained for several 
plate proportions and for one material. The results indicate plate 
strengths greater than those that have been found expert mentally on 
plates that do not satisfy straight-edge conditions. 


INTRODUCTION 


The determination of the load-carrying capacity of a plate sub- 
jected to loads in its plane depends upon a knowledge of the behavior 
of the plate in the postbuckled range. Postbuckling analyses of plates 
have for the most part been based on purely elastic considerations. 
However, the relatively thick plate elements used in modern aircraft 
structures may generally be expected to undergo plastic deformations 
prior to failure of the ccanponents that they constitute. Conseqiisntly, 
the theoretical determination of the loads that such plates can support 
requires the Incorporation of plasticity theory into a large-deflection 
postbuckling analysis. 

l»fe,ny authors have investigated the elastic postbuckling behavior 
of flat plates in conpressionj the more widely known of these investi- 
gations are references 1 to 11. The basic differential equations for 
a plate element undergoing large deflections are derived by Von Karmen 
in reference 1. In reference 2, Von Kerman introduces the concept of 
the effective width of postbuckled plates. Various approximate solu- 
tions for postbuckling behavior are presented in references 5^ 5^ 

and 6 by Cox, Timoshenko, Mar-guerre and Trefftz, and ftorguerre. 
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respectively^ where analyses are carried out hy energy methods. In 
reference 'J, Kromm and Marguerre obtain very accurate results at moder- 
ately exceeded buckling loads for simply supported, infinitely long 
plates in compression by extending the investigations of references 5 
and 6. An essentially exact solution for square plates in compression 
is presented by levy in reference 8, where the large-deflection equa- 
tions of reference 1 are solved to a high degree of approximation by 
means of Fourier series. In reference 9^ Koiter improves the results 
of reference 7 "to make them applicable far beyond buckling; in addition, 
results are presented for clamped plates. The effects of initial 
deviations from flatness for square plates are investigated by Hu, 
Lundqulst, and Batdorf in reference 10 and by Goan in reference 11 by 
means of the method of solution advanced in reference 8. In reference 10, 
the side edges of the plate are constrained to remain strai^t, whereas, 
in reference 11, the side edges are free to distort in the plane of the 
plate. 

The large niimber of investigations of the elastic postbuckling 
problem and the fact that solutions can be obtained only by approximate 
methods indicate to some extent the unwieldiness of the nonlinear, large - 
deflection equations involved in the postbuckling problem. Ihe lack of 
solutions to the corresponding inelastic problem appears to be a conse- 
quence not only of its large -deflect! on aspects but also of nonlinearity 
in the stress-strain relations. The inelastic problem is further conqjli- 
cated by the fact that some question exists as to the appropriate plastic 
stress-strain law applicable to the states of combined stress involved 
in the problem. 

In recent years, much work has been done in fonnulating variational 
principles for theories of plasticity. At the same time, great advances 
have been made in the development of hl^-speed ccanputlng machines. The 
use of a variational principle in conjunction with calculations performed 
on a high-speed computer appears to offer a feasible approach to the 
solution of the inelastic postbuckling problem. In the present paper, 
this approach is used to determine the behavior of an aluminum-alloy, 
sinply supported, square flat plate with edges constrained to remain 
strai^t, compressed beyond the (elastic) buckling stress into the 
plastic range. The variational principle used applies to the simple 
deformation theory of plasticity and the calculations were made on the 
Standards' Eastern Automatic Computer (SEAC) of the National Bureau of 
Standards . 

The results of the analysis are presented in the form of curves 
showing the variation of deformations, stress distribution, average 
con 5 )res 8 ive stress, and effective width with the applied imit shortening. 
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SYMBOIS 


E 


U 

V 


■b 


t) 


e 


e 


e 


cr 


f 


h 


p, q 
t 

tf 

U, V, V 


X, y, z 
^xy ' 


^xy 


ir 


^xy 


Yoiong's modulus for plate material 
secant modulus 

bending moments per unit length 
twisting moment per unit length 
strain energy of plate 
volume of plate 

plate dimension in x- and y-dlrection 
effective width of plate 

unit shortening applied to plate 
compressive buckling strain or shortening 

xmit displacement of plate side edges 

overall thickness of two-element plate measured 
between center lines of faces ^ 

integers 

thickness of homogeneous or solid plate 
thickness of face of t^ro-element plate 

displacement of point on middle smrface of plate in 
X-, y-y and z-direction, respectively 

plate coordinates (see fig. l) 

middle -surface shear strain 

twisting strain 

total shear strain in xy-plane 
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^eff 

P It II 


I 

n 

'^cy 

°cr 


o'eff 

Oy 

Oy 



Subscripts : 

t 

b 


effective strain. 





^x^y 


+ 


'xy 


middle-surface strains in x- and y-direction, 
respectively 

bending strains in x- and y-direction, 
respectively 

total components of strain in x- and y-direction, 
respectively 

nondlmensional lateral coordinate, 2 x/b 
nondimensional axial coordinate, 2 y/b 
Poisson's ratio for plate material 
average compressive stress 
conrpressive yield stress of material 


compressive buckling stress 

effective stress, + Oy^ - P^'^y ^^xy^ 

components of stress in x- and y-direction, 
respectively 

local average axial and lateral stresses 
local average shear stress 
shear stress in xy-plane 


top face of two-element plate 
bottom face of two-element plate 
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THEORY 

Statement of the Problem and Basic Assumptions 


The problem under consideration is to determine the behavior of a 
square flat plate ccm^ressed unidirectionally beyond the elastic buckling 
stress into the range where plastic yielding of the material takes place. 
All edges of the plate are assumed to remain in the original plane of 
the plate and to have vanishing bending moments. The plate is considered 
to be subjected to a uniform shortening by means of a pair of rigid, 
frictlonlesB loading platens; thus, the loaded edges remain straight and 
free of tangential stresses. The side edges - those parallel to the 
direction of loading - are assimied to be held straight but are free to 
translate laterally and are devoid of tangential stresses. The plate 
is further assumed (fig. l) to consist of two stress-carrying faces only, 
the stresses being constant throu^ the thickness of each face. The 
effects of finite transverse shear stiffness, however, are neglected. 

The choice of aspect ratio and edge conditions for the plate is 
based upon the observed behavior of an interior bay of a multiple-bay 
stiffened panel immediately after buckling. In general, nearly square 
buckles form in the skin, and the number of buckles in the longitudinal 
direction tends to persist in the postbuckllng range; thus, the buckles 
tend to remain nearly square. The behavior is somewhat different in 
the case of the tr^lly infinitely long plates analyzed in references 7 
and 9 vhere it is shown that the buckles tend to shorten in the direc- 
tion of loading as the load is increased. The restriction to square 
buckles made in the present paper is therefore considered to apply to 
panels that are long enough for nearly square buckles to be formed but 
not so long that the number of buckles can readily change. 

The two-element configuration assumed for the plate is incorporated 
into the problem in order to simplify the analysis by making it inde- 
pendent of the effects of plasticity throu^ the thickness of the plate. 

The question of what is the correct relation between stresses and 
strains imder combined stresses in the plastic region is at present 
quite unsettled. The literature contains analyses of plastic buckling 
(see refs. 12 and I5) that have been carried out, with evident success, 
on the basis of the simple secant -modulus deformation theory of plas- 
ticity. On the other hand, serious doubts exist concerning the legit- 
imacy of applying this theory when, as in buckling, the stresses may 
deviate sharply from so-called "proportional loading" - loading in which 
the conp)onents of stress at any point remain proportional to one another. 
A considerable controversy on this subject which will not be delved into 
in this paper has arisen. The course followed in this report is to use 
deformation theory but to consider plates that buckle elastically . In 
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such cases, where the posthuckled configuration has had an opportunity 
to hegin to develop elastically , the stress history in the plastic range 
■^^Quid not involve widespread, sudden deviations from proportional loading 
consequently, simple deformation theory may he used with some degree of 
confidence . 

In the present paper, Poisson's ratio is assumed to he equal to l /2 
in hoth the elastic and plastic domains) hence, material compressihility 
is neglected. According to the secant-modulus theory for an incOTpress- 
ihle material, the relationship between the instantaneous states of 
stress and strain at any point in the plate is given hy 


or 




1 

2 

1 

2 



‘xy 



s^xy 




(la) 



(lb) 


The quantity Eg is defined hy the relationship 

Eg 

^eff 

where Ugpf is an effective local stress given hy 


( 2 ) 


^eff = \f°x^ ^ ^y^ " '"x°'y + ^'^xy 

and €eff is ail effective local strain given hy 




£x *^v ^x^y + k 


( 5 ) 

(i^-) 
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The effective stress and strain are related hy the uniaxial stress- 
strain curve of the material. Bius, Eg can he determined hy entering 

the uniaxial stress -strain curve with either or 

Equations (l) are strictly applicable at any point in the plate 
only If the effective stress and strain at that point have never 
decreased while In the plastic range. The use of equations (l) for 
any stress history corresponds to a nonlinear elastic material whose 
uniaxial . loading and unloading stress-strain cixrves are identical. Hie 
present solution will be actually canried out for this hypothetical 
material by assuming equation (l) to apply always, and the applicability 
to the plastic material will then be assessed a posteriori by examining 
whether the effective strain does indeed increase monotonlcally once it 
enters the nonlinear range. 


Basic Equations 


In the present paper, the displacements u, v, and w of the 
middle surface (see fig. l) will be considered' as the basic unknoims 
to be sought. Since, after buckling, the displacement w can no 
longer be considered small in con 5 )arison with the plate thickness, the 
postbuckling problem requires application of finite-deflection theory 
to describe the relationships between strains and displacements. 


From reference 4, the expressions relating middle -surface strains 
and displacements are 


_ 8u , 

l/^w\ 

8x 

2\8x/ 


i/8w\ 

Sy 

2Wy 

= ^ 4 - ^ + 

oy 

ox 


2 

2 






8w 8w 
8x 8y 


(5a) 


The bending and twisting strains acting througjiout the thickness of a 
plate element are obtained from reference 4 as 




8^ 

8x^ 

8^w 




e = z 

y 8y2 

■y " = 

8x 8y 




(5t) 
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where %, in general^ is the disteince from the middle surface. In the 
case of the two-element plate, z is taken to be the distances 

from the middle surface to the median surfaces of the faces. 


With the use of the subscripts t and b to denote the top and 
bottom faces, respectively, the face strains can now be written as 


^t,b 


^ + 1 / 

'^w\^ + E 8^ 

8 x 2' 

'dx/ ” 2 

^ + 1/ 

'^wN^ 1 h 8^ 

8 y 2 * 

^ 8 y/ “ 2 ^y 2 


= ^ ^ ^ ^ + h 


^t,b 8y 8x 5 x 8y 


Sx 8y 
J 


(6) 


The stress-strain relations (eqs. (la)), in conjunction with the strain- 
displacement relations (eqs. (6)), give the stresses in the top and 
bottom faces as 


■^t,b 




H,b 




^t,b 


= i E. 


i 8 u , 

1 ^ 

+ i 

/ 8 w\^ 


2 Sy 

2 

\8xy 


1 ^ 

4. i 

/ 8 w\^ 


2 8x 

2 


+ 


8w 8w + 

\8y 

8x 

8x 8y 


1/^f 

2\Bx/ 


8 x Sy 


+ h /S^ , 1 8 £w\ 

+ h/ d^ , 1 S£w\ 
-^Sy 2 


>(7) 


Here Ea and E„ are considered to be functions of 


7 . 


^eff. 


— — S. g 

' Vt % yt 




l/T^ 


(8a) 


and 


^eff 


b 



2 2 




Tb 


•b yb 


(8b) 


respectively. 
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The local hending moments produced in the plate by the face stresses 
are simply 


n 4- 


-(o^t 

*V - -(oyt - “Vb)! 




^ V^t “ Vt,)2 


( 9 ^ 


J 


and the average stresses at any point in the plate are 


= 




•h 

■yh 


^ “ K'Vt * 

’'jiy “ i(^*yt * ’^^b) 




( 10 ) 


J 


The conditions of equilibrium in tiie x-, y-, and z-direction 
yield (see ref. 4 ) 


xy 


St. 

Sy 


= 0 


8t. 


xy _ 


8x 


= 0 




8x^ 


^ 8^ , - 8^ 


8^ 


8x 8y 8y2 8x2 ^ 8y2 ^ 8x 8y 


(11a) 

(lib) 

(11c) 


where 2t^ , the stress -carrying thickness of the two-element plate, is 

used in place of the thickness t of the homogeneous plate considered 
in reference 4 . 

The boundary conditions stipulated for the present problem are as 
follows: The prescribed unit shortening e requires that 
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The condition that the side edges remain straight is 


V 



= Constant 


The requirement that the side edges he free to translate is 




dx = 


0 


(15) 


(14) 


The condition that the tangential stresses vanish on each face is 
written 



Finally, the simple-support conditions stipulate that 
and 

Now, the differential equations (11), together with the boundary 
conditions (eqs. (12) to (17)), may he considered to constitute the 
complete statement of the problem in terms of the displacements u, v, 
and w. That is, through the use of equations (6) to (ll) and by spec- 
ification of the uniaxial stress-strain curve, the differential equa- 
tions (11) could, in principle, be reduced to a set of three differ- 
ential equations in u, v, and w. Similarly, the boundary conditions 
depend, implicitly, only on the displacements. An alternative formu- 
lation of the problem, less natural but more attractive, is presented 
in the next section by means of a variational principle. 
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Variational Principle 


The strain energy of stretching and bending for the tv?-o-eleinent 
plate is 


U = t^: 


^b/2 nb/2 / ne 

-b/2'-'-b/2 Wo 




<iev, + 






^7 


^t 




.b/2 

W2 / P 

'-h/2 

Wb/2 W 0 


de„ + 

^b ^b 


"^b 


de,,_ + 


0 


% 




0 


■^xy^ 


dx dy 


(18) 


It can be readily verified from equations (l) to (i^■) that 

°x ‘^^x ”^xy °^eff *^^eff ^^9) 

Thus, the strain-energy expression may be written 


U = t. 


>b/2 ob/2 

-b/2'^-b/2 




^efft 



dx dy 


( 20 ) 


As a result, in each face of the plate, the strain-energy density at a 
given point is simply the area under the uniaxial stress-strain curve 
up to the effective strain level at the point. According to the 
principle of mini mum potential energy for deformation theories (ref. 14- ) , 
a solution of the present problem renders the strain energy U a 
relative minimum. Stated precisely, 6U = 0 with respect to atlml sstble 
variations in u, v, and w. Py admissible variations are meant vsnia- 
tions that do not violate the geometrical boundary conditions of the 
problem^ in the present problem, the geometrical boundaiy conditions 
are given by equations (12), ( 15 ), and ( 16 ) . It is important to note 
that the variations in displacements need not satisfy the remaining 
so-called natural boundary conditions (eqs. (l4-), ( 15 ) , and ( 17 )). 

Althou^ the minimiim-potential-energy principle of reference l4- is 
formulated for small-deflection theory, its use in the present problem 
is valid. It is shown in appendix A, by means of the calculus of varia- 
tions, that minimization of the strain energy U with respect to 
admissible variations in displacements leads to satisfaction of the 
correct differential equations of equilibrium and the natural boundary 
conditions . 
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ffethod of Solution 

The solution of the present problem is found approximately by 
using the Rayleigh-Ritz procedure in conjunction with the variational 
principle discussed in the preceding section. Expressions for the 
unlcnowns u, v, and w are assiomed to be ,in the following form: 


00 00 

w = ^ ^ ^ 

^ • • • n~" 5 • • • 

00 CO 

u -ex + 

p=l q=0 
00 00 

V = f y + b^I y-nq cos 
p=0 q=l 


cos MX cos ^ (21) 

"b *b 

^ COB ^ (aa) 

t b 

^ sto SSS (25) 

b b 


These expressions satisfy all the gecmetrical boundary conditions of 
equations (12), (l5)^ and (16); e is the magnitude of the prescribed 
unit shortening and the unknown coefficients ^q^ ''^pq 

remain to be determined from the variational principle. In practice, 
the solution was limited to the determination of only the six unknown 
coefficients wt t , U]_q, ujq, f, Vqj_, and v-^ ^ . 

In the special case of elastic behavior, where is a linear 

function of egff and Eg = E everywhere in the plate, a solution for 

these coefficients by analytical minimization of equation (20) is 
feasible and is given in appendix B. In the plastic solution, however, 
it is necessary to introduce an appropriate uniaxial stress-strain 
curve in order to evalmte equation (20) and effect its minimization 
with respect to the six unknown coefficients. Such a minimization by 
analytical methods appears to be of prohibitive difficulty. Conse- 
quently, the course followed was to utilize a hi^-speed ccmiputing 
machine (SEAC) to carry out the required minimization process numeri- 
cally. Essentially, the procedure required repeated numerical evalua- 
tion of the strain-energy integral (eq. (20)) for systematically varied 
sets of the unknown coefficients. An exposition of the computational 
scheme is contained in appendix C. 


RESULTS AM) DISCUSSION 


The compressive stress-strain curve shown in figure 2, typical of 
24S-T aluminum alloy,, was used as the basic relationship between 
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o^ff and egff in the present investigation. Tlie postbuckling 

analysis was carried out for four different plate proportions, so 
chosen that the corresponding values of buckling stress were as indi- 
cated on figure 2; each plate was taken to have its initial buckling 
stress in the elastic range. 


As shown in appendix B, the elastic buckling stress for the two- 
element plate is given by 


■'cr 


5 vhy 


when Poisson's ratio equals 1/2, and the critical unit shortening is 



The postbuckling state for each plate depends only on ^/e^^, the 

ratio of the applied unit shortening to the critical strain. The final 
results for loads carried by two-element plates, with Poisson's ratio 
equal to l/2, are considered to apply approximately to solid plates 
having the same values of Ucj. and e/e^^. Ihe proportions of the 

eqiii valent solid plates are determined by the •usual formula: 


cTcr = 


E 


-r-Y 

5(1 - ^2) 


Displacements and Stresses 

By means of the numerical minimization process described in 
appendix C, approximations to the true values of the unknown coeffi- 
cients W]_ 2 .j ^10-> ^11^ v-| 1 in the expressions (21) 

to ( 25 ) for the displacements w, u, and v were determined for each 
of the four plate proportions considered and for -various values of the 
applied unit shortening ratio For a given unit shortening, the 

minimization process involved calculation of the strain-energy expres- 
sion ( 20 ) for systematically varied sets of the displacement coeffi- 
cients, so chosen as to ca-use the energy to decrease in magnitude 
contin^lally. For a given unit shor-bening, estimates of the true dis- 
placement coefficien-ts are then provided by those -values -that yield the 
lowest energy -value foxmd. With the use of numerical minimization, the 
accuracy of such estima-tes is limited by the fact that, in the neighbor- 
hood of its minimum, the energy is relatively insensitive to changes in 
the displacement coefficients. Nevertheless, a reasonably consistent 
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variation of the final coefficients with was found and is shown 

hy the faired c\irves of figures 3 an<i 4 for the two extreme plate geom- 
etries, plates 1 and 4, respectively. Also shown in these figures, for 
comparison, are the coefficients given hy the elastic solution of 
appendix B, which considers the same number of unknown coefficients. 
These elastic coefficients constituted the initial values used in the 
iterative process for minimizing the energy in the plastic range. As 
can he seen from figures 5 14-, the largest differences between the 

elastic and plastic displacements occur for the in-plane displace- 
ments u and v; the deflection coefficient w-| is very nearly the 
same for hoth the elastic and plastic cases . 

From the values of displacements, stresses can he found hy use of 
equations (7) • The lateral distrlhution of axial stress 9^^, as deter- 
mined hy equation (lO), is shown in figure 5(a) for three different 
cross sections of plate 1 at a value of e/e^^^ of 5 ; similar distribu- 
tions for plate 4, at a value of e/e^j, of 6 , are also shown in fig- 
ure 5(’b). These results for the two plates actually correspond to the 
same unit shortening, since = ^(®cr) 2 _- Figures 5(c) and 5 (d) 

show the distributions of average lateral stress Oy at several longi- 
tudinal cross sections of each of the two plates. It may he of interest 
to note that, in contrast to the results of figure 5 ^ the corresponding 
elastic solution of appendix B yields axial and lateral stress distribu- 
tions (eqs. (bio)) that are independent of the x- and y-coordinate, 
respectively. 

In the elastic solution, conditions of equilibrium in the x- and 
y-direction are satisfied exactly at each point of the plate. The 
extent to which in-plane equilibrium is satisfied hy the plastic solu- 
tion may he measinred in a gross sense hy checking the closeness of the 
values of the resultant axial force at various lateral cross sections; 
also, the lateral stress distributions may he examined to see whether 
they produce essentially zero resiiltant lateral force. The results of 
such equilibrium checks are illustrated in figure 5 ^ where the average 
stress values are indicated for each cross-section distribution. It is 
seen that, although the average lateral stress is close to zero in all 
cases, there remains some discrepancy between the magnitudes of the 
average axial stress at the various stations. Based on the mean value 
of average axial. stress for the three cross sections in each case, the 
percentage spread between the largest and smallest values of average 
axial stress amounts to 4.8 percent for plate 1 and 7-^ percent for 
plate 4. It may be of interest, however, to ccmpare these discrepancies 
ivLth those that would be obtained by calculating plastic stresses on the 
basis of the elastic displacements that were lised in the first cycle of 
iteration in the minimization process. The average axial-stress valiies 
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obtained on this basis are indicated in figure 5 "by the values enclosed 
in parentheses. As can be seen, the spread between the largest and 
smallest values of average axial stress for each of the plates is sub- 
stantial, 32.4 percent for plate 1 and 29-5 percent for plate 4. Thus, 
a significant improvement was achieved by the minimization process. 

The strain distributions of plates 1 and 4 and their variations 
with s/e^j. were subjected to a detailed examination in order to deter- 
mine whether plastic unloading ever occiarred. It was found that, up to 
the largest value of e/e^j. considered in the present study, no such 

unloading did in fact occur. That is, the magnitude of the effective 
strain £ 0 ff in each face continually increased with increasing e at 

all points in the plate where plastic yielding occurred) at some points 
(for example, the midpoint of the loaded edge) ^eff lioes decrease 

immediately after buckling but this unloading occurs in the elastic 
range. 


Relation Between Average Applied Stress and Unit Shortening 

The load-carrying capaciiy and the stiffness of a plate after buck- 
ling are determined by the relationship between the average compressive 
stress aa_-^r and the applied unit shortening. The load-carrying capacity 
of the plate at a given shortening is singjiy the product of the average 
compressive stress and the cross-sectional area whereas the stiffness of 
the plate after buckling is related to the slope of the curve of average 
stress against unit shortening. 

As has been shown, the resultant axial force obtained in the present 
solution varies somewhat, depending on the location of the cross section. 
An appropriate unique value of the average applied stress may be most 
conveniently determined by means of the following energy considerations . 

The strain energy (eq. (20)), which was minimized for various values 
of e/scj., must equal the total external work. Therefore, 



whence 


■^av 


de Y 


(24) 


where the plate volume is V = 2tj.b^. Since the minimization process 
described in appendix C yields values of U/EV directly, the use of 
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relation ( 24 -) to determine quite convenient. (The same 

relationship was used In the elastic analysis of reference 9 >) 

Figure 6 shows the variation of U/EV with e/e^P for each of 

the four plates investigated. The required differentiation of these 
curves was effected hy means of the graphical procedure described in 
reference I5, and the resultant curves of against e are shown 

in figure 7 "the four plates. For ccraparison, the curves determined 
by the elastic solution of appendix B are shown for each case. 

It should be emphasized that these results are certainly approxi- 
mate by virtue, among other things, of the fact that only a limited 
nunber of coefficients were used in the expressions for the displace- 
ments . An estimate of the extent of the inaccuracy may be made by 
con5)aring the essentially exact relationship between average stress and 
unit shortening found by levy for the elastic case with that given by 
the approximate elastic solution of appendix B, which considers the 
same nuimber of coefficients as were irsed in the plastic solution. This 
comparison is made in figure 8 , which shows cfay/hcr function of 

e^e^j, for these two solutions. Thus, it is reasonable to expect that the 
curves of figure 7 reliable only up to values of ^/^cr about 

this limiting value is noted by tick marks on each of the four curves . 

A reasonable procedure for correcting the curves of figure 7 would 
be to reduce them at each val\ie of e/e^.j. by the corresponding ratio 

of the two curves of figure 8 . This reduction has been made in figure 9 } 
the curves constitute the final estimates of the present paper for the 
relationship between average compressive stress and unit shortening for 
the four plates studied. 

An Interesting feature of the results of figure 9 is that the 

average compressive stress carried by each plate does not have a maxi- 

mum value in the range of unit shortening considered. However, the lack 
of definite mathematical maximums in the average compressive-stress — 
unit-shortening curves of figure 9 should, perhaps, not be overempha- 
sized in view of the flatness of these curves and the approximate nature 
of the present solution. More significant is the fact that the curves 
indicate load-carrying capacities substantially greater than those that 
have been found experimentally for plates that do not conform to strai^t- 
edge boundary conditions. Such test results (unpublished) have been 
obtained by R. A. Anderson and M. S. Anderson of the Langley laboratory 
for plates supported laterally in a V-grooved fixture similar to that 
used in the classical tests by Schuman and Back (ref. I6) . Plates 
supported by V-grooved fixtures are not, of course, constrained to have 
strai^t (in -the -plane) side edges after buckling, and, in addition, 
the V-grooved fixture does not entirely prevent out-of -plane displacements 
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of the side edges once initial buckling takes place. Needham, in 
reference IT, gives the test results found hy a number of investigators 
for the maximum strength of con^ressed square tubes. Again, the plate 
elements of square tubes do not satisfy the boundary conditions of 
straight side edges beyond buckl in g, and, as is remarked by Needham, 
failure occurs with passage of the buckle pattern through the corners. 

A comparison between these test resiilts and the results of the present 
paper can be made by taking from figure 9 ^ nominal maximum average 
stress eqml to the stress at the large unit shortening of 0.01 

(the highest unit shortening for which ccsnputations were carried out) . 
Figure 10 shows the relationship between 0max/°’cr ‘^cr/'^cy 

in this fashion from figure 9 also gives the corresponding c\irve 
determined by test results on a variety of materials, including 
24S-T aluminum, for square tubes (ref. I 7 ) and for plates in V-grooved 
fixtures tested by Anderson and Anderson. (ihe results for both types 
of tests fall essentially on a single curve.) 

It appears, then, that maintaining straight-edge boundary conditions 
leads to hi^er load-cari*ying capacities, as indicated by the results of 
the present analysis. Hiis conclusions is bolstered by the experimental 
results found by Botman in reference 18 , which are also shown in fig- 
ure 10. Botman, continuing the Investigations started by Besseling 
(ref. 19 )^ used a jig that divided a wide plate into three strips by 
means of a series of opposed knife edges running longitudinally. This 
type of fixture represents more closely the boundary conditions of the 
present investigation than do the V-grooved fixtures or the square tubes 
and, as can be seen frcm. figure 10, leads to higher plate strengths. Of 
course, the outside edges of the outer bays still do not conform to the 
straight-edge condition^ one may conjecture, then, that still hl^er 
strengths would be achieved from tests on plates with additional bays. 


Effective Width 

The so-called "effective width" of plate can be readily calculated 
frcm the curves of figure 9 by means of the relationship 





where Ogg is the stress obtained from the stress-strain curve for the 
value of e that yields a given value of variation of 

bg/b with e/ecr is shown in figure 11 and is compared with the effec- 
tive width calciilated on the basis of Levy's "exact" elastic solution 
for the square plate. Ihls comparison shows that, at the higher unit 
shortenings, the effective widths become somewhat greater than those 
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given by elastic theory. The early dip of the present results helov 
the elastic curve is a consequence of the fact that plastic strains 
occur in the interior of the plate at unit shortenings that are still 
in the elastic range. 


CONCLUDING REMARKS 


A theoretical analysis has been made of the behavior of a simply 
supported, square flat plate with strai^t edges con^iressed beyond the 
buckling load into the plastic range. Approximate solutions were carried 
out numerically for one material and for four different plate proportions; 
the curves for average compressive stress against unit shortening thus 
obtained were corrected in a rational fashion to account for the limited 
number of degrees of freedcxa used in the analysis. 

The method of analysis, involving the application of a variational 
principle in conjunction with a high-speed computing machine, may be 
applied to plates with different boundary and loading conditions and of 
different material properties. The number of degrees of freedom assumed 
in the displacement functions may be increased, subject to the capacity 
of the computing machine. The chief difficulty that was encountered in 
the actual numerical solution was the slowness of satisfactory convergence 
to minimum energy values. The development of Improved methods of effecting 
numerical minimization of nonlinear functions of many variables would be 
hi^ly desirable for future application of the method of analysis used in 
this investigation. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 25, 195^* 
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APPENDIX A 

VERIFICATION OF VARIATIONAL PRINCIPLE 


The purpose of tiiis appendix is to sLow that the variational 
condition 


SU = 0 (Al) 

where U is given hy aqiiation (20) , is valid for finding a solution to 
the present prohlem. The variation of U is taken with respect to 
admissihle variations in u, v, and w, that is, those that do not vio- 
late the geometrical houndary conditions (12), (13) > and ( 16 ) . It will 
he shorn that the differential eqmtions (U) and the natural boundary 
conditions (eqs. (l4), ( 15 ), and (I 7 )) are consequences of equation (Al) 

If the strain-energy density 


'^eff 


'0 


cTeff <i€eff 


is denoted hy > equation (20) may he written 

^h/2 p-h/2 

U = tf / / (F^ + Pb) dx dy 

'■^-h/2 ^-h/2 


where the subscripts t and h refer to the top and bottom faces, 
respectively. Then, the variation of U with respect to u, for 
example, is 


oh/2 ph/2 

= t^ / / ^^u^t ^u^h^ (-^) 

^ -h/2^ -h/2 


Similar expressions hold for variations with respect to v and w. 
Now, 


6^F 


= ^ 


Su^x + 


8f 


&u^ + 


8f 

87 . 


xy 


^u7xy 


(A3) 
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andj again, similar expressions hold for S^F and 6^. Since 
<1F = Cgfp d£g£.^, it follows from equation (I9) that 


8 f 


5 f 

8f 

87 




= a 


X 


= cr„ 


= T 


xy 


xy 


From the strain-displacement relations (6), it follows that 

8 


8u^x = 


(A 4 ) 


Hence 


8u€y = 0 


Vxy = 8u(^“) = 


iJJ . tj. 


ph/2 

rh/2 r 

^-h/2 

’-■h /2 1 






dx dy 


After integration hy parts and with the use of expressions (lO), the 
variation becomes 


5u^U = 2tf 


Ph/2 

x=h/2 

Ph /2 

/ ^ 6u 

dy + 

^-h/2 ^ 

"^-h/2 

yL=-\>f 2 


y=8/2 

y=-h/2 


dx - 


_ dT^\ 

■h/2 '^-h/S \ Sx 8y / 


5 u dx dy 
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= +:k 

"2 


On the other 


But, 5u = 0 at X = ±— (Bince a = xe— at x 

^ 2 V 2 

hand, 5u is arbitrary at y = hence, the condition 5^U = 0 
requires satisfaction of the boundary conditions 

= ° 

as -well as of the equilibrium, equation (lla). 


8 a 


8t 


0 


8x 8y 

Similarly, calciilation of 5.^U leads to 


SyU = 2t^ 


~ob/2 

y^b/2 

pb/2 

/ a- 8v 

'^-b/2 ^ 

dx + 

^-b/2 

y=-b/2 ^ 


x=b/2 

x=-b/2 


dy - 


Pb/2 

Pb/2 / 

^ -b/2 ^ 

-b/2 \ 


6vd.dy 


Hence, setting S^U equal to zero reqiiires satisfaction of equation (Ub) 


8ct 8t 

_Z + _J 2 L = 0 

8y 8x 
- 


Also, since 6v is arbitrary at x = ±^, it follows that 


T^(±|,yj = 0 


On the other hand, by virtue of the straight-edge condition (eq. (I 3 )), 


6v must be constant at y = ±^. Hence, 


^b/2 

-b/2 




:,±|) = 0 


since the constant value of 6v at each side edge may be arbitrary. 

Thus far, it has been shown that the equilibrium equations (lla) 
and (Ub), as well as the natural boundary conditions (lij-) and (15)., 
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are consequences of condition (Al) . It will now be shown that the 
remaining equilibrium equation (lie) and natural boundary conditions (I7) 
will follow from the condition 6^U = 0. 

From equations (6 ) , 

^ ^ h(bv) , h h^(bv) 

^ ^t,b 5 x - 2 Sx 2 

§ g = ^ 5 ( 6 w) + h 8^(Sw) 

^ yt,b 8y 8y 2 ^y2 

5 -y = ^ 8(Sw) ^ 8(5w) + 8^(Sv) 

^ ^t,b Sx Sy 8y 8x 8x by 


Hence , 

6..U = t. 


-'w'' 


'^b/2 nb/2 
-b/2 ^-b/2 


(“ 


+ a 


\ 8w 8(&w) 
^/8x 8x 




Vb) 


8w 8(5 w) ^ ^ 8(5w) 
8 x 8y Sy Sx 


(^x. 


h 5^(Sw) 

t - Sx2 


(oy - Oy. ^ + (t^ - T \h >■ dx dy 

\^t yby2 5y2 \^t ^b/ 8x by 


Through the \ise of relations (9) and (lO) and after integration by parts, 
the variation becomes 


5wU 


■''b/2 

-b/2 






Sy 3x 


8v 


'-■b/2 ^ Sy 


y=b /2 

y=-b/2 

x=b/2 

x=-b/2 


dx - 


3y Sx 
-b/2 


Sy / 

dv 


x=b/2 


^b/2 


x=-b/2 
x=b/2 


ay + 


Pb/2 

-b/2 


I * * 


r 


b/2 


x=-b/2 
y=b/2 

'-b/2 0*|y=-b/2 


dy - 


-^b/2 


MyS 






dx - 


-b/2 " 5y 

pb/2 pb/2 


y=b/2 

y=-b/2 


dx + 


'-b/2 ^ -b/2 




2tf 0 ^ || + 4t, V ^ ^ . 2 ^ dx dy 
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Now, if 6,^11 is to vanish for all admissible values of 
and 5 then it is readily seen that satisfaction of 

oy 

requires that 


5w_, 


6 


3w 


= 0 




0 


and 




- 2 


8x 8y 8y2 V Sx^ 


- 8^ , 


2t. 


5^w 


^y2 " Sx Sy^ 


Therefore, it has now been shown that all the natural boundary conditions 
(eqs. (1^1-), ( 15 ) , and ( 17 )) and the differential equations of equllibrixmi 
(eqs. (lla), (lib), and (lie)) follow from the condition that the first 
variation of the potential energy of the plate must vanish for all 
admissible variations in the displacements u, v, and w. 
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APPENDIX B 
ELASTIC SOLUTION 


The elastic posthuckllng behavior of the two-element plate can be 
obtained analytically by applying the Rayleigh-Ritz method to the 
elastic counterpart of the variational problem formulated in the present 
paper . 


The strain energy of stretching and bending of the plate is given 
by equation (I9) except that it is now understood that the relationships 
between the components of stress and strain are equations (la) with the 
local modulus Eg replaced by Young's modulus E. Substitution of these 
elastic stress -strain relations into equation (18) and integration to the 
final strain state yield 


U = ^ Et. 


^b/2 ob/2 
-b/2 '^-b/2 


2 2 ’'^t^ 


I Xt) y-b x-b y-b 


XYb 


dx dy 


(Bl) 


Since the displacements u^ v, and w can be related to the strains 
through the Inrge-def lection strain-displacement relations (6), the strain 
energy can be written as a function of the unit shortening e and the 
undetermined coefficients appearing in the displacement expressions (2l), 
(22) and (25). As mentioned earlier, both the plastic and the elastic 
solutions of the present paper are limited to the determination of only 
six of the unknown coefficients appeaning in the displacement functions. 
That is, the displacements are assumed to be 


u 


V 


-ex + B^UpQ + U-] 1 cos 
fy + b^VQp + v^p cos 


2zti^gin irtx 
b / b 

^>in ^ 

b / b 






(B2) 


w = bw-i T cos — cos ^ 
b b 


J 


where e is the applied unit shortening and the remaining coefficients 
are to be determined by application of the Rayleigh-Ritz method. The 
assumed expression for w is exact at initial buckling for a square 
plate and is assimed to be reasonably accurate in the early postbuckling 
range . 
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Through the .use of relations ( 6 ) and (B 2 ) in eq.uatlon (B 1 ) , the 
strain-energy density can he integrated over the stress-carrying area 
of the plate to give 


U = i Etf-h^ <j e^ - ef + f^ + I + f) + ^ 


( 2 rtUio)^ + ( 2 rtvoi)‘ 


^(2:iniio + 2:rvoi)w3j_2 + — 


( 2 j£Uii )2 + ( 2 jtvii )2 


2^(2jniii + 2 «vii)vii 2 + |(2itun) (2jtvn) + 


(B5) 


The conditions for obtaining the values of the coefficients that 
minimize the strain energy for a given value of the unit shortening are 

M = 8 u ^ 8 U ^ Bu ^ 5 U ^ 8 U ^ Q 
Bf 8 u]_o ^■'^01 ^^11 ^^11 


After the operations indicated hy the minimizing conditions are performed, 
the following six simultaneous eqmtions are obtained for evaluating the 
six unknown coefficients: 

- - + f + ^ = 0 

2 16 ^ 

^2 p 

2 ( 2 rtuio) - ^ wpp"^ = 0 
o 

2 ( 2 jtVQi) - ^ wjj_^ = 0 

f (2ituii) + ^(2itvpi) - ^1^2 

4 A 4 

^ cr 2 

^( 24 Uii) + ^( 2 jtvii) - ^ wii^ 


= 0 


= 0 


(b4) 


W' 


11 


:l(_e + f) - i 
2 2 


l/ 2 rtU 3 _o + 2 itVQp 


+ 2 rtu-| 1 + 2 jtv-| 1 + 







0 
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Solution of the first five of equations (b 4) yields 

'ir 


f = ^ - 2 


16 


2«Uqo = 2jtVQ2_ - 


w 2 

16 


2 

2rtUii = 2rtVii = 2^ WijL^ 


CB5) 


J 


Substitution of relations (B5) into the last of equations (Bh) gives 


W' 


11 


-e + 


4 ^th\2 «2 2 

5(t; 


= 0 


(b6) 


In general, W 2 _q ^ 0; hence, 

5Uy 


-e + 


£ 

4 




(BT) 


Now, when W]_q vanishes, the criterion for initial buckling is estab- 
lished; that is, the critical strain or shortening becomes 


= |(f )‘ 


(b8) 


and, in turn, the compressive buckling stress is simply 


a,, . it EfskV 
3 Vb y 


As a consequence of equation (b 8) , equation (BT) yields the nondlmensional 
center deflection of the plate as 


wii' 


2 _ ^ 


‘+ / 


'cr 


(B9) 


and, hence, the -six coefficients are completely determined for any value 
of the plate unit shortening and cross-section geometry. 

By virtue of the results (B5) and (B9 ) , the displacement expres- 
sions (B2) can be substituted into the stress-displacement relations ( 7 ) 
(with Eg = E) to obtain the stresses present in the faces of the two- 
element plate. In terms of average stress (see relations (lO)), the 
results of this substitution yield 
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It should he mentioned that, although the stress expressions (BIO) were 
obtained on the basis of two-element-plate geometry and a Poisson's ratio 
of 1/2, a rederivation for the case of a solid plate, with an arbitrary 
Poisson's ratio, yields the identical results for the stresses. 

It is of Interest to note that the approximate elastic solution 
(eq.s. (bio)) for the stresses turns out to satisfy the eq.ullibri\m equa- 
tions (lla) and (Ub) exactly at each point of the plate; the third 
equilibrium equation (eq. (lie)) is, of course, not satisfied. A measure 
of the accuracy of the approximate elastic solution is afforded by 
figure 8 , which compares the approximate relationship found between 
%v/°”cr ®/®cr with the exact relationship found by Levy. 

Equations (B 5 ) and (B 9 ) determined in this appendix are used as 
initial values in the minimizing procedure for determining the solution 
to the plastic-plate problem. (See appendix C.) 
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APPENDIX C 


PLASTIC SOLOTION 


The variational problem is conveniently formulated in nondimensional 
terms as foUovs: Let | = 2x/b and = 2y/b; and let 

= I 

Since the four quadrants of the plate behave identically^ the energy 
integral (20) may be modified into the nondimensional form 



Here is given by equations (8) in -terms of and 

■which in turn are found from equations ( 6 ), (21), (22), and (25) to 
depend only on the spatial coordinates | and 7^, on the applied unit 
shortening e, and on the unknown nondimensional parameters 

2itUio 2jtVQi 

2Trui3_ 2rtv^ 

^ ^ ^11 

The ratio Sll fixes the magnitude of the critical shortening e„„ and 
b 

hence specifies the plate geometry. The set of six nondimensional 
parameters (C2) are to be so determined as to minimize F for given 

values of — and e. 
b 

Since analytical minimization of F is not feasible when nonlinear 
stress-strain relations are involved, recourse is had to a numerical 
minimization process in conj-unction -with -the use of a hl^-speed computing 
machine (SEAC) . The minimization process was effected by means of, 
essentially, the so-called "method of steepest descents." (See, for 
example, refs. 20 and 21.) 
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The basic idea of the procedure may be described in general terms 
as follows: Consider a function F(^xp,X 2 , • • -x^) . The set of 

n independent parameters may be conveniently denoted by the n-component 
vector XjL (i = 1, 2, 5^ • • • n) . Ihe value of Xj_ that minimizes F 

is sought. An initial trial vector Xj_(^) is assumed^ and the gradient 

of F, that is. — , is calculated at Xj(^). The direction - is 

8xi' ^ Sxi 

then the direction of steepest descent of the function F; the func- 


tion F 


8xp 


is then evaluated for various positive values 


of 5 in an effort to find the value 6 that minimizes F 




When this value of S is found (presianably, approximately ) , a new direc- 
tion of steepest descent is determined by evaluation of the gradient of 

F at the point - 6 xp(*^M . The process is continued 

Sxj^ L J 

until satisfactory convergence is obtained to the lowest possible value 
of F. In the present problem, the six parameters in the set (eq.. (C2)) 
play the role of the con 5 >onents of the vector Xj_. 


In the application of the method of steepest descents, the basic 
procedure outlined above was modified in certain respects. First of all, 
each evaluation of F was performed approximately by means of numerical 
integration. EJach face of the quarter plate was divided into a grid of 
100 sqmres, the integrand in eqiiation (Cl) was computed at each grid 
point, and the integral was evaluated by applying Simpson's rule twice - 
once for the integration in the |-direction and once for the T|-direction. 
For this calculation, the function G was represented, piecewise, by 
polynomials detennlned from the stress-strain relationship of figure 2. 

The evaluation of the gradient of F was performed on the basis of 
a finite-difference approximation to each of the six partial derivatives 

required. Thus, the determination of the six components of — neces- 

OXi 

sitated seven evaluations of F - one at the particular set of starting 
values of the independent variables and one each for a small increment 
in one of the six variables. 

In any given cycle, the value of 6 in the starting point Xj_ - 5 

8xj_ 

of the succeeding cycle was found by passing a parabola throu^ the three 
points F(o), F(6) , and F(26) given by the values of F at x^, 

x^ - 6 and Xj - 25 . Here 6 was chosen to be of some conven- 

^ 8xi ^ 8xi 

lent magnitude, preferably of the order expected for 6. 
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It is pointed out in reference 20 that, in practical application, 
the method of steepest descents tends to furnish successive approxima- 
tions to the minimizing vector that zigzag tovard the true minimum 
rather than approach it in a smooth fashion. Consequently, the following 
procedure was introduced in an attempt to speed up convergence: After 

two successive cycles of minimization in the direction of the negative 
gradient, a third cycle of a different nature was inserted. In this 
extra cycle, minimization was performed not in the direction of the 
negative gradient, hut rather in the direction determined hy the differ- 
ences between the last-ohtained approximations to the unknowns and the 
approximations of two cycles before. Ihus, a round of three successive 
approximations to the minimum vector proceeded as follows: 

(1) Given x-?^: find ^ minimize 

ZhCj_ V / 

find obtain x^^^^ = 

(2) Repeat step (l), starting with this procedure gives 


n 



to 


(5) Let 




^i 


n+2 


= ^ 


n+1 


- 6 


'nfl 


AF 

AXj_ 


Sj_ = x^^ - x^^^; minimize F 
Xi^5 = - E^2^±' 



to obtain 


(4) Repeat steps (l), (2 ), and (5) as many times as are necessary 
to obtain satisfactory convergence to the minimizing value of x^^. 


It may be remarked that in practice the quantities and were 

Axj_ 

normalized; that is, each component of a set or Sj was divided 

AXj 


by the absolute value of the largest component. 


It is clear that a tremendous amount of numerical calculation is 
Involved in the application of this procedure to the present problem; 
only very hi^ cociputing speeds make such a proced'ure feasible. 


For each value of and e, the initial values assumed in the 

b 

iteration process were taken from the elastic solution of appendix B. 
Convergence to the minim\mi energy level was considered satisfactory if 
several successive cycles yielded the same value of energy to within 
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four signlficaxit figures. About thirteen cycles of iteration, involving 
approximate ly 35 minutes of con^juting time for the SEAC, wefe generally 
required for convergence. In several instances, particularly at the 
highest unit shortenings, convergence was found to he iiiqjracticahly slow; 
in these few cases, the minimum energy level was estimated hy evaluating 
it for values of the displacement coefficients extrapolated from curves 
such as are shown in figures 5 4. The final plots for energy against 

unit shortening obtained for the foiir plates are shown in figure 6. 

It is of interest that the total difference in energy between the 
first and last cycles of iteration for eny one case was always very small 
(of the order of 1 percent) . Consequently, as far as the energy is 
concerned, a fair estimate could be found on the basis of the elastic 
displacement coefficients. However, as is discussed in the body of the 
present paper, a reasonably accurate estimate of the stress distribution 
requires application of the minimization process. 
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(a) Axial stress 


Plate 1; (Td) Axial stress 


Plate 4> 


. ksi 0.0272 .0224 

' ''av 


,0336 .00875 


(c) Lateral stress, 3y. Plate 1; (d) Lateral stress, 5y. Plate k} 

- 5 . -S- = 6. 

cr ®cr 

Figiore 5 .- Axial and lateral stress distributions for first quadrant of 
plates 1 and 4. (Values in parentheses refer to average stress based 
on elastic displacements.) 











cvj rO 5^ 





Figure 8.- Comparison of approximate and exact elastic solutions for the 
average compressive stress carried by a buckled plate. 
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Figure 



. - Variation of average compressive stress with unit shortening 

(adjusted) . 
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O Computed values of 


at one-percent 


unit shortening from present theory 
□ Square-tube tests (ref. 17) 

O V- groove tests (unpublished) 

A Knife-edge tests, 3- bay plates (ref. 18) 


Figure 10 .- Coii5)arison of experimen'fcal and calculated reBults for com 

pressive strength of flat plates. 




M>01 - SS-W-r - i»iau»l-VOVN 



,4 L L 

.3 ,4 



J I I 

8 .9 1.0 


Figitre U 


■1.A.U n ui_Lkwi 

results of Levy's 


^ nea/^ *f*-KVTm Tvr*Ckd=i'n+ c/^T Ti+-f r^i-\rrrrM3 -tj-t +Vi 

WW S. W .*. I ^ ,1 «-< s*. ^ 4. W4JJL WI>^W4i^ V W V.^ ?V Ul-L 

exact elastic analysis (ref. 8). 


liACA TN 5368 



